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Measurement of the Two-Jet Differential Cross 
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Abstract 

A measurement is presented of the two-jet differential cross section, d 3 a/dE T dr]idr] 2 , 
at center of mass energy y/s = 1800 GeV in pp collisions. The results are based on an 
integrated luminosity of 86 ptT 1 collected during 1994-1995 by the CDF collaboration 
at the Fermilab Tevatron collider. The differential cross section is measured as a 
function of the transverse energy, E T , of a jet in the pseudorapidity region 0.1 < |?7i| < 
0.7 for four different pseudorapidity bins of a second jet restricted to 0.1 < \r] 2 \ < 3.0. 
The results are compared with next-to-leading order QCD calculations determined 
using the CTEQ4 and MRST sets of parton distribution functions. None of the sets 
examined in this analysis provides a good description of the data. 
PACS numbers: 13.85.Rm, 12.38. Qk 
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Jet production in proton-antiproton collisions results predominantly from hard 
interactions between two initial state partons. Theoretical developments in both per- 
turbative next-to-leading order (NLO) and parton shower Monte Carlo calculations 
permit calculation of many QCD jet processes with theoretical uncertainties small 
enough to allow detailed comparison with measured distributions fl]. In this paper, 
we present a measurement of the dijet differential cross section that provides more 
precise information about the initial state partons than has been probed by previous 
CDF measurements of inclusive jet transverse energy 0, total transverse energy [BJ, 
and dijet mass All previous measurements showed an excess of events at high jet 
energies when compared to the QCD prediction based on standard sets of parton dis- 
tribution functions (PDFs). One explanation for this excess is a larger than expected 
number of high momentum partons, particularly gluons, in the proton |6j . While 
those measurements provide cross sections averaged over a wide range in their vari- 
able, in this analysis we reduce the region over which averages are taken by measuring 
the cross section for four separate ranges. This provides more detailed information 
about the cross section shape. Previous measurements of the dijet differential cross 
section have been performed by the CDF and D0 |8j collaborations with smaller 
data samples. The present measurement places new constraints on the parton distri- 
butions of the proton. 

Jet production rates are usually expressed in terms of the transverse energy, Et-i 
and pseudorapidity, r), of the jets, where rj is related to the polar angle 9 relative to 
the proton beam line by rj = -ln[tan(0/2)]. At leading order in QCD, the proton, 
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p, and anti-proton, p, momentum fractions, x\ and x 2 , carried by the two colliding 
partons can be expressed as 

x x = %{e^ + e" 2 ), x 2 = %{e~^ + e~ m ). (1) 

V s V s 

Here 771 and r] 2 are the pseudorapidities of the two jets, yfs is the center of mass 
energy of the colliding hadrons and Ej- is the transverse energy of the leading jet. 
For a fixed Et and 771, one can probe higher x values by selecting events in which the 
second jet has a larger 7/ 2 value. For a given x we have four measurements at what 
are effectively different values of Q 2 the square of the four-momentum transfered in 
the interaction, calculated by 

Q 2 = 2£!coshV(l -tanhr/*), v * = -(rji - r) 2 ) . (2) 

The four distributions in this analysis allow us to measure the cross section on a 
surface in the x-Q 2 phase space whose shape is sensitive to the predictions of different 
PDFs. 

The constraint on the parton distributions at high x comes mainly from prompt 



photon production in pp or pA collisions from WA70 and the E706 [1C] experiments 
and inclusive jet data from the Tevatron @. The data do not constrain the parton 
distributions very well at high x. The higher statistics of this measurement together 
with the multiple cross section measurements at different Q 2 for approximately the 
same x provide a precise set of data which can be used to determine improved sets 
of PDFs. The current measurement, based on data of an integrated luminosity of 
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86 pb 1 from 1.8 TeV pp collisions taken during the 1994-1995 Fermilab Tevatron 
collider run, covers the range 0.05 ~ X\ ~ 0.8. 



The CDF detector is described in detail in [nj] . In this analysis we utilize the cen- 
tral, plug, and forward calorimeters. The central calorimeter covers the pseudorapid- 
ity range < 1.1. It is segmented into projective towers of size At] x A0 = 0.1 x 0.26, 
where <fi is the azimuthal angle in radians. The plug (1.1 < \r}\ < 2.4) and forward 
(2.4 < \t]\ < 4.2) calorimeters are segmented by approximately 5° in and 0.1 in r\. 
The event vertex is resolved to within 1 mm along the z axis, using time projection 
chambers surrounding the beam pipe. 



A cone algorithm with cone radius R = J (A0) 2 + (A77) 2 = 0.7 is used to identify 



jets | 12| . Transverse energy is defined as Et = Esm9, where E is the scalar sum of 
energy deposited in the calorimeter towers within the cone and 9 is the angle formed 
by the event vertex, the beam direction, and the cone center. Our data sample 
consists of events collected by on-line identification of at least one jet with transverse 
energy above trigger thresholds of 20, 50, 70, and 100 GeV at integrated luminosities 
of 0.091, 2.2, 11, and 86 pb -1 , respectively. The bin widths in Et were chosen to be 
larger than the measurement resolution on Et and to ensure sufficient statistics in 
the bins. 

In this analysis we use events with at least two jets of Et > 10 GeV of uncorrected 
energy. We consider events in which the -Sr-weighted centroid of at least one of the 
two highest Et jets is in the range 0.1 < \r]\ < 0.7. This "leading" jet is required 
to deposit more than 40 GeV Et, prior to corrections, in the central calorimeter. 
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In addition, the centroid of the second leading jet is required to be in the region 
0.1 < \r)\ < 3.0, and the primary event vertex must be located within ±60 cm of the 
nominal interaction point. Poorly measured events and background from cosmic rays, 
beam halo, and detector noise are removed by requiring that total energy recorded 
by the detector be less than 2000 GeV and E T /y/EE^ < 6V~GeV, where E T is the 
missing transverse energy and T,Et is the scalar sum of the total transverse energy. 

In this analysis, we evaluate the E T spectrum of the leading jet for the following 
four r] bins of the second leading jet in the event: 

0.1 < 1 772 1 < 0.7, 0.7 < 1 772 1 < 1.4, 

1.4 < 1 772 1 < 2.1, 2.1 < 1 772 1 < 3.0. 

The 772 ranges were chosen to place regions of reduced response (due to gaps between 
detectors) within single bins while at the same time maintaining a sufficient number of 
events in the bins. Both jets are included in the distribution for the 0.1 < |t7 2 | < 0.7 
bin if each satisfies the requirement 0.1 < |?7| < 0.7 and Et > 40 GeV. 

Since the calorimetric response varies as a function of 77, we determine the trig- 
ger response separately for each 772 bin. The trigger efficiency was measured using 
overlapping E T regions for the different trigger thresholds. For the 20 GeV trigger 
threshold, for which no lower E T trigger was available, the second jet in the event was 
used to determine the trigger efficiency. For the four trigger thresholds, the trigger 
efficiency was found to be greater than 90% for jets of E T greater than 40, 82, 105, 
and 130 GeV. 

The measured jet Et must be corrected for calorimeter non-linearity and loss of 

8 



energy in the gaps between calorimeters. In addition, the measured jet Et spectrum 
must be corrected for the smearing effect caused by the resolution in the measured 
jet E T . We simultaneously correct all these effects with the procedure used in our 
previous measurement of the inclusive jet Et spectrum ||. For the central rj bin 
(0.1 < 1 771 1 < 0.7) at 40 GeV, the correction to the measured Et is approximately 
4%, while the correction to the measured cross section is about 19%. The correction 
to the cross section increases to 70% for the bin 2.1 < |ryi| < 3.0. The corrected cross 
section values are given in Tables [1] and and plotted in Figure [I]. 

The systematic error on the measurement of the jet cross section is dominated 
by the uncertainty in the measurement of the jet Et magnified by the steep slope 
of the Et spectrum. Although the same sources of uncertainty contribute to the 
cross section of each E T bin, the uncertainty depends on the local slope of the E T 
spectrum. The systematic uncertainties were evaluated as in References |§ and [|i~3 



The uncertainties include: charged hadron response at high pt (h pt); calorimeter 
response to \ow-p T hadrons (1 pt); ± 1 % on the jet energy of the absolute calibration 
of the calorimeter (esc); jet fragmentation functions used in the simulation (frag); 
± 30% on the underlying event energy in the jet cone (uevt); detector response to 
electrons and photons (e/ph); and modeling of the detector jet energy resolution 
(cres). The resolution on the measured r/ causes events to migrate between adjacent 
bins. In the highest rj bin, the gap between the plug and forward calorimeters results 
in decreased r] resolution and has the effect that more events migrate out of the 
bin than into it. To compensate for this effect, we have applied an E^-dependent 
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Table 1: The measured dijet differential cross sections for 0.1 < 1 772 1 < 0.7 and 
0.7 < 1 772 1 < 1-4. The differential cross section is given for the average E T of the bin. 
The statistical and systematic errors are shown as a percentage of the central value. 
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Table 2: The measured dijet differential cross sections for 1.4 < |r] 2 | < 2.1 and 
2.1 < 1 772 1 < 3.0. The differential cross section is given for the average E T of the bin. 
The statistical and systematic errors are shown as a percentage of the central value. 



11 



correction which is less than 8% in all bins. The effect was studied by breaking it 
into two components, the resolution on the measured 77 (77 res) and a systematic shift 
in the reconstructed rj (rj sh). It is included in the systematic error by looking at 
the result on the cross section when doubling and halving the correction. Bins for 
which events were collected using triggers with uncorrected energy greater than 20 
GeV (J20), 50 GeV (J50) and 70 GeV (J70) were assigned 4, 2 and 2 percent errors 
respectively, associated with prescaling. An overall luminosity uncertainty (norm) of 
4 percent is added in quadrature with these. The sources of systematic errors are 
listed in Tables [3] through |6| as percentages of the central Et value for each Et and 
77 bin. In general the percent error increases as 772 increases. 

In Figure |2], the difference between the fully corrected two-jet differential cross 
section and the predicted cross section is divided by the predicted cross section and 
plotted as a function of the leading jet E T for the four 77 ranges of the second jet. 
The theory predictions were calculated using the NLO calculation of the JETRAD 



program [14] with the PDFs indicated. The calculations use a renormalization scale 
\x = E™ ax /2 with Rs e p = 1.3, where R sep is a measure of the maximum separation 
between the cones of two jets that are merged into one. The error bars represent 
the statistical errors, while the shaded bands represent one standard deviation of the 
systematic error, which is correlated for all the different Et values. The data are 
compared to the predicted cross section obtained using the PDF set CTEQ4M [RJ. 
The solid curve shows the expected results when using CTEQ4H J |J , and the dashed 
curves show the results when using the PDF set MRST 
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Table 3: The systematic errors for the 0.1 < 7] 2 < 0.7 bin given as a percentage of 
the central value. The E T values are specified at the bin average. A reference to the 
sources of systematic errors is given in the text. 



13 



<E T > 


e/ph 


uevt 


frag 


esc 


cres 


1 pt 


h pt 


rj sh 


rj res 


norm 


J20 


J50 


J 70 


tot 


43.1 


2.6 


15.0 


8.3 


4.2 


6.1 


7.7 


2.7 


0.6 


0.3 


4.0 


4.0 


2.0 


2.0 


21.5 


49.5 


2.6 


12.3 


8.5 


4.1 


5.4 


7.6 


3.1 


0.5 


0.3 


4.0 


4.0 


2.0 


2.0 


19.6 


58.5 


2.7 


9.5 


8.6 


4.1 


4.6 


7.4 


3.6 


0.5 


0.3 


4.0 


4.0 


2.0 


2.0 


17.8 


75.2 


2.8 


6.6 


8.5 


4.1 


3.4 


7.0 


4.6 


0.5 


0.3 


4.0 


4.0 


2.0 


2.0 


16.2 


94.0 


3.0 


5.0 


8.4 


4.3 


2.6 


6.5 


5.7 


0.4 


0.2 


4.0 


4.0 


2.0 


2.0 


15.7 


106.2 


3.1 


4.5 


8.3 


4.5 


2.3 


6.2 


6.4 


0.5 


0.2 


4.0 


0.0 


2.0 


2.0 


15.1 


119.1 


3 2 


4.1 


8 2 


4 8 


2.2 


5 Q 


7 2 


5 


2 


4 





2 


2 


15.2 


131.9 


3.3 


3.9 


8.1 


5.1 


2.1 


5.6 


7.9 


0.5 


0.2 


4.0 


0.0 


0.0 


2.0 


15.4 


150.0 


3.5 


3.8 


8.1 


5.5 


2.3 


5.2 


9.0 


0.6 


0.3 


4.0 


0.0 


0.0 


2.0 


16.0 


173.2 


3.8 


3.7 


8.2 


6.1 


2.6 


4.9 


10.3 


0.7 


0.3 


4.0 


0.0 


0.0 


0.0 


16.9 


206.9 


4.2 


3.6 


8.7 


7.1 


3.3 


4.7 


12.2 


1.0 


0.5 


4.0 


0.0 


0.0 


0.0 


18.9 


260.5 


5.1 


3.5 


10.4 


8.9 


4.5 


5.0 


15.2 


1.7 


1.1 


4.0 


0.0 


0.0 


0.0 


22.9 


313.7 


6.0 


3.4 


13.2 


10.9 


5.4 


6.0 


18.2 


2.6 


2.3 


4.0 


0.0 


0.0 


0.0 


27.7 


373.9 


7.1 


3.2 


17.9 


13.5 


6.7 


8.1 


21.7 


3.9 


4.7 


4.0 


0.0 


0.0 


0.0 


34.6 



Table 4: The systematic errors for the 0.7 < r] 2 < 1.4 bin given as a percentage of 
the central value. The E T values are specified at the bin average. A reference to the 
sources of systematic errors is given in the text. 
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Table 5: The systematic errors for the 1.4 < i]2 < 2.1 bin given as a percentage. The 
Et values are specified at the bin average. A reference to the sources of systematic 
errors is given in the text. 
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Table 6: The systematic errors for the 2.1 < r] 2 < 3.0 bin given as a percentage. The 
E T values are specified at the bin average. A reference to the sources of systematic 
errors is given in the text. 
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The observed excess of events at high Et values in the inclusive jet cross sec- 
tion measurement may be explained within the framework of conventional QCD by 
exploiting the relatively weak restriction on the gluon density at high x §. The 
CTEQ4 PDFs use a more flexible parameterization of the gluon density at high x 
than is present in other sets. The CTEQ4 set of PDFs include the inclusive jet data 
from the Tevatron. The CTEQ4HJ PDF gives a higher weight to the inclusive jet 
data while still maintaining agreement with the other data sets used in the fit. 

The MRST set of PDFs is based in a wide range of deep inelastic scattering 
data and has an improved treatment of heavy flavors and prompt photon production 
than do previous MRST sets. The main constraint upon the gluon at high x comes 
from prompt photon production from the WA70 |9| and E706 ]IIJ data. The set 



MRST(gf) was derived assuming that there is no initial state partonic transverse 
momentum (< > = 0); this does not lead to a good fit for the prompt photon 
data from the E706 experiment. The set labelled MRST(gj) was derived by allowing 
non-zero < > while maintaining reasonable agreement with the WA70 data. The 
MRST(gj) set has <k^> = 0.64 GeV. These two sets represent the extreme values of 
<kr> that yield reasonable agreement with the data used in the fit. The set labelled 
MRST represents the preferred set from the global analysis and has < kx > = 0.4 GeV. 
The covariance matrix for the dijet cross section is 

Vij = 8ijOi(stat) + T<lLiVi(sysk)crj{sySk), 

where dij = 1(0) for % = j(i ^ j) ai(stat) is the statistical uncertainty in bin i and 
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Ci{sysk) is the systematic uncertainty, k, on bin i. The sum is over the 13 sources 
of systematic errors listed above and over all the bins in each of the four 77 bins. 
We calculate the \ 2 from \ 2 — ^ijAj(V r ~ 1 )j, ? -A J -, where A3 is the difference between 
the data and theoretical prediction for bin i. The average of the upper and lower 
errors is used when calculating the x 2 - The x 2 /dof values for different PDFs are 
presented in Table ^. Although the cross sections predicted by the MRST PDFs 
are lower than the data by 20%, they have similar \ 2 values to those predicted with 
CTEQ4M. This is because the systematic errors allow a correlated shift in the data 
which makes only a small contribution to the total \ 2 ■ Predictions whose shape 
matches that of a correlated systematic error will give reasonable \ 2 values provided 
that the normalization between the data and prediction are within a few standard 
deviations. The probability of describing the data with the PDFs used in this analysis 
is less than 1% in all cases. 
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2.88 
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4.49 







Table 7: The \ 2 1 dof between the data and the prediction of Giele, et a!., for different 
PDFs. The fit to the data has 51 degrees of freedom. 

In summary, we have measured the differential cross section for dijet production 
in pp collisions with one jet restricted to the pseudorapidity region 0.1 < \rji\ < 0.7 for 
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four different pseudorapidity bins of a second jet restricted within 0.1 < 1 772 1 < 3.0. By 
allowing the pseudorapidity of the second jet to vary through 0.1 < \rj\ < 3.0, we are 
able to map out the cross section over the available kinematic phase space and provide 
a differential cross section that more tightly constrains the parton distributions of the 
proton than in measurements previously reported by us. The measurement provides 
more precise information about the parton distributions of the proton in the high x 
region, an area which is not well constrained, and will provide useful input to QCD 
global fits. The resulting improved sets of PDFs will help to further enhance our 
knowledge of the structure functions of the proton. 
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Figure 1: The differential cross section for dijet events as a function of transverse 
energy, Et, and pseudorapidity, rj, of one jet, for 4 ranges in the pseudorapidity of 
the other jet. The results are compared with QCD predictions using different parton 
distribution functions. 
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Figure 2: The differential cross section for dijet events as a function of transverse 
energy, E T , and pseudorapidity, r], of the leading jet, for four ranges in the pseudo- 
rapidity of the second leading jet. The results are compared with QCD predictions 
using different parton distribution functions. The statistical error is represented by 
the error bars while the correlated systematic error is shown as the shaded band. 
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